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The  purpose  of  this  study  was  to  examine  the  response  of  deeply-burled 
tunnels  to  the  effects  of  high  stress  dynamic  loading.  In  the  Phase  I 
portion  of  the  study,  attention  was  concentrated  on  the  effects  of  geologic 
discontinuities  (e.g..  Joints)  on  tunnels,  with  particular  emphasis  on  the 
Interactions  between  joints  and  backpacking.  In  Phase  II,  emphasis  was  on 
the  response  of  Integral  tunnels  in  hard  and  soft  rock  geologies,  and  on  the 
implementation  of  a coupled  finite  element/finite  difference  method  (called 
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INTRODUCTION  AND  SUMMARY  OF  RESULTS  OF  PHASE  I 


1 . 

All  rock  masses  contain  large  numbers  of  nominally  planar 

discontinuities,  such  as  joints,  faults,  cracks,  and  bedding 

planes.  These  surfaces  usually  have  properties  which  are 

quite  different  from  the  basic  medium,  and  there  is  both  ex- 
1 2 

perimental  and  theoretical  evidence  to  suggest  that  the 
effects  of  these  discontinuities  should  be  one  of  the  import- 
ant considerations  involved  in  designing  survivable  under- 
ground structures. 

In  extreme  cases,  the  sliding  of  material  along  joints 
or  other  discontinuities  could  cut  off  a buried  structure 
(by  guillotining),  or  it  could  split  a tunnel  lengthwise,  as 
shown  in  the  sketches  below.  In  less  extreme  cases,  the  sit- 
uation shown  in  Figure  1 may  develop,  where  portions  of  the 
tunnel  wall  attempt  to  project  into  the  tunnel.  The  tunnel 
liner  must  resist  such  motion.  Backpacking,  of  course,  makes 
it  easier  for  a tunnel  liner  to  accommodate  small  displace- 
ments without  damage. 


In  the  Phase  I study,  interactions  of  the  type  shown  in 
Figure  1 between  deep ly- bur icd  tunnels  and  in  situ  joints  in 
granite  when  loaded  by  1 kb  stress  waves  were  investigated 
using  both  analytical  models  and  two-dimensional  finite 


Loading  Wave 


difference  numerical  code  calculations.  The  study  consisted 
of  (1)  the  development  and  use  of  analytical  solutions  to 
simplified  models  to  predict  if,  and  where  slippage  of  joints 
occurs  near  a tunnel,  (2)  the  development  of  a model  which 
provides  a conservative  estimate  of  the  ability  of  a back- 
packing system  to  resist  and  accommodate  block  motions,  and 
(3)  the  generation  of  seven  exploratory  numerical  solutions 
(using  the  Lagrangian  WAVE-L  finite  difference  code)  to  exa- 
mine the  t ime- resolved  wave  propagation  interactions  between 
various  lined  and  unlined  tunnels,  with  and  without  joints 
present  in  the  surrounding  rock.  The  conditions  for  these 
numerical  solutions  in  Phase  I are  listed  in  Table  1.  All  of 
the  solutions  were  two-dimensional,  plane  strain  cases  treat- 
ing a 1 kb  planar  step  wave  incident  on  an  infinite  cylindri- 
cal tunnel.  Models  were  developed  for  use  in  the 
numerical  solutions  to  describe  the  hydrodynamic,  elastic, 
and  failure  properties  of  the  backpacking  material.  Also, 
constitutive  properties  were  formulated  for  dilatent  rock 
joints. 

Two  important,  related  conclusions  can  be  drawn  from  the 
results  of  the  Phase  I study; 

1 . Significant  block  motions  and  joint  slippage  in  the 
immediate  vicinity  of  a tunnel  will  be  primarily  activated  by  ■ 

stress  relief  waves  emanating  from  the  opening.  This  activa-  : 

tion  of  joint  motions  is  most  likely  to  occur  near  the  tunnel  ) 

springlines,  but  it  involves  only  blocks  which  are  near  the  4 

tunnel.  There  may  be  small  amounts  of  joint  slippage  in  the  | 

free  field  before  the  stress  wave  reaches  the  tunnel,  but  the  I 

tendency  for  joint  slippage  is  greatly  increased  near  the  tun-  I 

nel  due  to  local  stress  relief.  (It  is  of  course  evident  that 
large  relative  displacements  can  occur  across  major  fault  j 
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Table  1.  Conditions  for  Numerical  Solutions 
in  Phase  I Portion  of  Study 


Case 

Medium 

Joints 

Backpacking 

Liner 

A 

Granite 

None 

None 

Uniined 

B 

Granite 

~3'  spacing, 
frictionless , 

_L  to  wave  front 

None 

Unlined 

C 

Granite 

None 

5'  simulated 
backpacking 

1"  rigid 
steel 

D 

Granite  w/no 
tensile  failure 

3'  spacing,  diag. 
to  wave,  dilatent 
w/friction 

None 

Uniined 

E 

Granite 

3'  spacing,  diag. 
to  wave,  dilatent 
w/friction 

None 

Unlined 

F 

Granite 

3'  spacing,  diag. 
to  wave,  dilatent 
w/friction 

None 

1"  rigid 
steel 

G 

Granite 

3'  spacing,  diag. 
to  wave,  dilatent 
w/frictlon 

5'  cellular 

concrete , 

P -1500  psl 
c ^ 

3/8"  rigid 
steel 

1 kb  planar  step  wave  loading 
applied  to  upper  boundary  of  grid 

t t { T t V' 

7 


60  ft 


1'  nominal  zoning 
used  in  all  problems. 
Tunnels  were  at  4000' 
depth. 
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planes  or  long  planar  surfaces  with  low  shear  resistance,  with 
or  without  tlie  presence  of  a tunnel,  but  we  were  not  concerned 
with  these  large-scale  phenomena  in  this  study.) 

2 . Backpacking  and  friction  are  very  effective  in  limit- 
ing these  local  block  motions  for  0.5  - 1.5  kb  stress  loads  on 
deeply-buried  tunnels.  Figures  2 and  3 show  estimates  (obtain- 
ed from  an  analytical  model)  of  the  thicknesses  of  backpacking 
required  to  stop  a block  of  length  L which  has  been  accelerat- 
ed into  a tunnel  by  0.5,  l.C,  and  1.5  kb  stress  loads.  The 
backpacking  resists  the  block  penetration  with  a crushing 
stress  of  1900  psi  in  uniaxial  strain.  After  the  backpacking 
has  sustained  9%  strain,  it  begins  to  lock  up,  and  stresses 
rise  sharply  above  1900  psi.  In  our  conservative  model,  we 
assumed  that  this  condition  would  result  in  tunnel  failure. 
Figure  2 shows  the  extreme  case  of  frictionless,  or  open 
joints.  Even  in  this  worst  case,  11  inches  of  backpacking  is 
sufficient  to  stop  a 10-ft  long  block  accelerated  by  a 1 kb 
wave.  25  inches  of  backpacking  are  required  for  a 1.5  kb  wave. 
Figure  3 illustrates  the  effectiveness  of  frictional  forces 
acting  along  the  joints  bounding  the  moving  block.  With  a 
nominal  frictional  shear  stress  of  100  bars  (which  roughly 
corresponds  to  a "wet"  joint),  only  4.5  in.  and  10  in.,  res- 
pectively, are  required  to  stop  a 10-ft  long  block  accelerated 
by  1 kb  and  1.5  kb  waves. 
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Backpacking  Thickness  Required  Before  Block  Motion  into  Tunnel  Stops  (inches) 


Backpacking  Thickness  Required  Before  Block  Motion  into  funnel  Stops  (inches) 


Figure  3.  Backpacking  Thickness  Required  to  Stop  a Block 
Accelerated  into  a Tunnel  by  Stress  Loadings  of 
0.5,  1.0,  and  1.5  kb  (100  bar  frictional  stress 
on  sliding  Hurf.icCH  of  block) 
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SUMMARY  OF  RESULTS  OF  PHASE  IT,  AND  RECOMMENDATIONS  FOR 

CONTINUING  EFFORT 


2 

2.1  PHASE  1 1 SUMMARY 

In  Phase  II,  primary  attention  was  given  to  liners 
in  softer  geologic  media,  as  exemplified  by  NTS  tuff,  and  to 
integral  liners  (i.e,,  liners  without  porous  backpacking). 

This  Phase  consisted  of  the  following  main  efforts: 

o Development  of  a model  and  associated  numerical  tech- 
nique which  allow  shear  cracks  to  explicitly  form  and 
propagate  in  the  Lagrangian  finite  difference  code, 

o Exploratory  numerical  solutions  to  examine  effects  of 
dynamic  stress  loading  on  lined  and  unlined  tunnels  in 
tuff,  and  to  examine  the  effects  of  explicit  shear 
crack  failure  (as  contrasted  to  conventional  plastici- 
ty failure)  on  the  tunnel  response, 

o Development  of  the  LAP  code,  a linked  finite  differ- 
ence/finite element  method,  for  coupled  analysis  of 
the  response  of  thin  liners  in  geologic  materials,  and 

o Coupled  numerical  solutions,  using  LAP,  to  analyze 
the  response  of  steel  liners  (with  and  without  con- 
crete backing)  in  tuff. 

The  conditions  for  the  six  numerical  solutions  performed 
in  Phase  II  arc  summarized  in  Table  2.  As  in  Phase  I,  all  of 
these  solutions  were  two-dimensional,  plane  strain  cases  in 
which  the  tunnel  was  subjected  to  a 1 kb  planar  step  wave. 

Some  of  these  cases  were  exercising  or  demonstrating  new  as- 
pects of  the  analytical  methods  for  treating  buried  structures 
for  the  first  time  (i.e.,  the  explicit  shear  crack  failure 
model  and  the  coupled  LAP  code) . 

The  results  of  the  Plmsc  II  efforts  are  described  in 
the  balance  of  this  report.  The  most  important  of  the 
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Table  2.  Conditions  for  Numerical  Solutions 


in 

Current  Phase  II 

Portion  of 

Study 

Case 

Medium 

Fracture  Model 

Liner : 

Thickness  Material 

H 

Granite 

Explicit  shear 
cracks 

Unlined 

J 

NTS  Tuff 

Plasticity 

Unlined 

K 

NTS  Tuff 

Plasticity 

1-ft 

Concrete 

L 

NTS  Tuff 

Explicit  shear 
cracks 

1-ft 

Concrete 

M 

NTS  Tuff 

Explicit  shear 
cracks 

2-in 

Deformable 

steel 

N 

NTS  Tuff 

Explicit  shear 
cracks 

2-in 

+ 

2-ft 

Deformable 

steel 

Concrete 

1 kb  planar  step  wave  loading  applied 
to  upper  boundary  of  grid 


t * t < t t 


T 


60  ft 


1'  nominal  zoning 
used  in  all  problems. 


Remarks 


LAP  code 
(linked  FE/ 
FD)  sol 'ns 
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findings  are  the  following: 

1 .  Shear  cracks  (and  presumably  tensile  cracks  also) 
which  form  near  the  tunnel  as  a consequence  of  the  dynamic 
loading  can  lead  to  non-uniformities  in  velocities  near  lined 
tunnels,  and  hence  to  stress  concentrations  in  liners.  These 
stress  concentrations  will  need  to  be  taken  into  consideration 
in  the  design  of  integral  liners. 


2 .  The  LAP  code  provides  a practical  method  for  analyz- 
ing the  coupled  dynamic  response  of  thin  liners  and  the  sur- 
rounding medium.  The  finite  difference  portion  of  the  method 
is  used  to  analyze  free  field  motions  and  the  medium  response 
(including  joint  slippage  and  generation  of  explicit  cracks), 
while  the  finite  element  portion  is  used  to  analyze  the  res- 
ponse of  the  liner.  This  combination  avoids  the  necessity 
of  using  small  cells  (and  small  time  steps)  to  resolve  action 
in  a relatively  thin  liner  in  a purely  finite  difference 
analysis . 


3 .  Two-inch  steel  liners  in  10-ft  dia  tunnels  in  tuff 
subjected  to  1 kb  stress  wave  loading  will  experience  stresses 
which  are  large  enough  to  cause  inelastic  deformation.  Diame- 
tral strains  of  about  2-3  inches  develop  shortly  after  the 
tunnel  liner  is  engulfed  by  a 1 kb  step  wave.  Accelerations 
(relative  to  the  center  of  mass  of  the  tunnel)  associated  witn 
these  distortions  are  of  the  order  of  1000-3000  g's. 

2.2  RECOMMENDATIONS 


Using  the  coupled  LAP  method,  realistic  liner  systems,  as 
well  as  new  liner  concepts,  should  now  be  analyzed  to  predict 
their  response  in  various  media  and  with  various  design  para- 
meters. LAP  is  able  to  treat  fairly  complex  2-D  designs,  and 
to  generate  dynamic  response  predictions  which  would  be  very 
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useful  in  making  preliminary  evaluations  of  new  concepts,  in 
examining  dynamic  failure  mechanisms,  and  in  helping  to  opti- 
mize design  variables.  The  following  specific  work  is  re- 
commended : 

1.  New  design  concepts  - The  experimental  evaluation  of 
new  design  concepts  is  both  expensive  and  time  consuming. 
Small-scale  experiments  are  useful  for  preliminary  evaluations, 
but  there  are  uncertainties  and  limitations  in  these  tests 
relative  to  boundary  conditions,  material  property  simulation, 
and  in  observing  potentially  important  dynamic  phenomena  on 
such  small  scale.  LAP  code  solutions  offer  an  alternative  or 
supplementary  approach  for  preliminary  evaluations.  Such 
solutions  are  particularly  valuable  in  that  they  indicate  the 
failure  modes  or  weaknesses  of  new  designs,  thereby  suggesting 
avenues  for  improvements. 

2.  Dynamic  failure  mechanisms  - The  Phase  I results  in- 
dicate that  relative  displacements  across  joints  will  not 
usually  be  larger  than  a few  inches  under  conditions  of  inter- 
est. (We  refer  here  to  the  local  motions  which  are  discussed 
in  the  Phase  I summary  in  Section  1 above.)  Such  motions  can 
be  accommodated  by  backpacking,  but  they  become  a much  more 
serious  concern  with  integral  liners,  which  must  resist  any 
motion  which  would  lead  to  significant  relative  displacements 
near  the  tunnel.  In  resisting  such  motions,  large  local  con- 
centrations of  forces  applied  to  the  liner  may  develop,  as  in- 
dicated in  the  sketch  below: 


In  addition  to  force  concentrations  due  to  motions  along 
preexisting  joints  or  cracks,  the  work  in  Phase  II  with  expli- 
cit crack  propagation  suggests  that  cracking  which  is  induced 
by  the  stress  wave  loading  near  the  tunnel  may  also  produce 
non-uniformities  in  loading  on  the  liner.  The  dynamics  of 
such  force  concentrations  and/or  non-uniformities  need  to  be 
considered  in  liner  design. 

We  had  hoped,  during  Phase  II,  to  be  able  to  analyze  one 
or  more  cases  such  as  shown  by  the  sketch  on  the  preceding 
page.  However,  we  did  not  get  that  far  after  completion  of 
the  LAP  code  implementation.  In  addition,  we  felt  that  such 
analyses  should  be  performed  after  we  could  treat  both  explicit 
shear  cracking  and  tensile  cracking  simultaneously  in  the 
Lagrangian  code.  We  recommend  that  this  additional  modeling 
capability  be  added  now,  and  that  the  local  interactions  of 
displacement  discontinuities  with  tunnel  liners  then  be 
analyzed  to  determine  the  stress  concentrations  in  liners. 

3.  Parameter  optimization  - The  LAP  method  offers  a 
fairly  efficient  way  to  conduct  series  of  analyses  which  can 
be  used  to  guide  the  optimization  of  design  variables  so  as 
to  obtain  the  maximum  useable  tunnel  volume  for  a given  ex- 
cavated volume.  Survivability  criteria  can  be  specified,  for 
example,  in  terms  of  liner  fracture,  or  maximum  allowable 
distortion,  velocity,  or  acceleration.  Solutions  of  this 
type  can  be  obtained  at  modest  cost,  and  the  results  used  to 
construct  curves  showing  sensitivity  of  liner  response  to  the 
design  variable,  medium  properties,  etc. 
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MATERIAL  PROPERTIES  FOR  NUMERICAL  SOLUTIONS 


5.1  TUFF 

3.1.1  Equation  of  State 

Emphasis  'n  the  Phase  II  portion  of  this  study  has  been 
on  tunnels  in  relatively  soft  geologic  materials.  NTS  tuff 
was  chosen  as  being  representative  of  these  softer  materials. 
Fairly  extensive  constitutive  property  data  are  available  for 
this  material,  as  summarized  in  Ref.  3.  These  data  show  that 
the  properties  of  NTS  tuff  vary  considerably  with  depth,  wat- 
er content,  etc.  Because  of  this  variability,  and  because 
the  peak  overpressures  experienced  by  the  tuff  in  our  solu- 
tions would  be  only  about  1 kb,  a relatively  simple,  hydro- 
static equation  of  state  with  constant,  but  different,  loading 
and  unloading  moduli  was  adopted,  based  on  the  data  in  Ref.  3. 
We  believe  the  important  characteristics  of  deeply-buried  tuff 
in  the  stress  range  of  concern  here  are  adequately  represented 
by  this  approach. 


A sketch  of  the  behavior  of  the  hydrostatic  equation  of 
state  used  in  loading  and  unloading  is  shown  below: 
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As  indicated,  the  tuff  loads  with  a bulk  modulus,  = 90  kb. 

It  unloads  (and  reloads)  with  a bulk  modulus,  = 116  kb. 

This  behavior  was  chisen  to  model  hysteretic  compaction  of  the 
material.  Unloading  from  4 kb,  for  example,  results  in  a 
volume  reduction  of  1%,  which  is  consistent  with  experimental 
hydrostatic  loading  curves  and  residual  volume  compaction  data 
for  several  NTS  tuffs. 

A constant  shear  modulus  during  loading  and  unloading  was 
used,  based  on  experimental  bulk  modulus  and  longitudinal 
sound  speed  data  from  Ref.  3.  The  values  used  were  28  kb  for 
the  shear  modulus  and  = 8500  fps. 

3.1.2  Failure 

The  failure  properties  of  the  tuff  were  modeled  in  two 
ways:  The  first,  which  we  refer  to  as  the  plasticity  model, 

is  the  conventional  method  in  which  failure  is  assumed  to  mani- 
fest itself  as  a yielding  phenomenon  occurring  on  a plastic 
yield  surface.  The  rationale  for  this  model  is  that  the  mater- 
ial fails  by  developing  numerous  small  cracks  which  simulate 
plastic  behavior.  The  second  failure  model,  wliich  we  refer  to 
as  the  explicit  shear  crack  model,  represents  a refinement  of 
the  plasticity  approach.  This  method  allows  the  formation  of 
explicit  shear  cracks,  along  with  appropriate  stress  adjust- 
ments . 

Material  failure  data  for  NTS  tuff,  taken  from  Ref.  4, 
are  shown  in  Figure  4.  These  data  give  the  maximum  shear 
stress,  T,  which  any  plane  embedded  in  the  material  can  exper- 
ience without  failing,  as  a function  of  the  normal  stress,  a^, 
on  that  plane.  For  use  in  the  numerical  model,  the  data  in 
Figure  4 were  approximated  by  the  function 
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Figure  4.  Maximum  Shear  Stress  Which  Can  Be  Withstood 
Across  a Plane  vs  Normal  Stress  on  Plane, 
NST  Tuff  (from  Ref.  4) 


where 


a . 


T = .16  kb 

oo 

c = .115  kb 

a = .465  kb 
o 


Failure  Model 


The  yield  stress  for  the  conventional  plasticity  failure 
model  was  obtained  from  eq.  (1)  by  interpreting  the  normal 
stress,  as  the  pressure,  and  the  shear  stress  as  being 

proportional  to  the  yield  stress  at  that  pressure.  The  pro- 
portionality constant  was  chosen  to  be  2,  which  (under  con- 
ditions of  uniaxial  strain)  results  in  yielding  being  initiat- 
ed at  the  stress  level  where  shear  cracks  first  appear.  Since 
the  numerical  solutions  have  a uniaxial  strain  far  field 
boundary  condition,  a direct  comparison  can  thus  be  made  be- 
tween the  two  failure  models. 


)licit  Shear  Crack  Failure  Model 


In  this  model,  material  fails  by  developing  explicit 
shear  cracks  whenever  the  shear  stress  across  any  plane  ex- 
ceeds the  value  given  by  eq . (1).  After  failure,  appropriate 
stress  adjustments  are  made  in  the  cell  to  reflect  the  exist- 
ence of  the  cracks.  The  explicit  shear  crack  model  also  re- 
quires a post-failure  function  expressing  the  maximum  shear 
stress  across  an  already  existing  crack  as  a function  of  the 
normal  stress  across  the  crack.  In  the  absence  of  better  data, 
we  used  half  the  shear  stress  provided  by  eq.  (1)  for  this 
post-failure  condition.  Note  that  in  isotropic  materials, 
shear  cracks  form  in  pairs:  when  the  shear  stress  across  one 
plane  exceeds  that  required  for  a fracture  to  form,  there  is 
always  another  plane  across  which  this  criterion  will  also  be 
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satisfied.  Thus,  in  the  code  solutions,  no  cell  ever  develops 
a single  shear  crack;  they  always  come  in  pairs.  In  field 
plots  obtained  from  the  numerical  solutions,  the  orientation 
of  the  shear  cracks  is  indicated  by  line  segments  which  are 
drawn  in  each  failed  cell.  If  more  than  two  shear  cracks  are 
formed  in  any  cell,  the  plasticity  failure  model  is  thereafter 
adopted  for  the  material  in  that  cell. 

Explicit  tensile  cracking  has  also  been  programmed  in  the 
WAVE-L  Lagrangian  code  in  connection  with  another  application, 
but  the  combination  of  explicit  tensile  and  shear  cracking  is 
not  currently  allowed  in  the  code. 

3.2  GRANITE 

For  the  one  granite  solution  in  Phase  II  (Case  H) , a lin- 
ear elastic  material  model  with  bulk  modulus,  K = 512  kb, 
shear  modulus,  G = 293  kb,  and  density,  p = 2.7  g/cm^  was  as- 
sumed. The  shear  failure  surfaces  for  wet  and  dry  granite 
joints,  as  determined  in  the  Phase  I effort,  are  reproduced 
in  Figure  5.  For  Case  H,  we  assumed  that  the  shear  failure 
surface  for  intact  granite  is  the  same  as  for  a dry  granite 
joint . 

3 . 3 CONCRETE 

A simple,  linear  elastic  equation  of  state  was  used  to 
represent  the  concrete  liners,  with  a bulk  modulus,  K = 270  kb, 
and  a shear  modulus,  G = 138  kb.  The  concrete  was  not  allow- 
ed to  fail  in  these  solutions. 
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Figure  5.  Shear  Failure  Surfaces  for  Wet  and 
Dry  Joints  in  Granite 


4.  UNLINED  AND  CONCRETE- LINED  TUNNELS 


4.1  UNLINED  TUNNEL  IN  GRANITE 

Case  H involved  an  unlined  tunnel  in  granite.  The  granite 
was  modeled  as  being  linear  elastic,  failing  only  when  shear 
cracks  develop.  Tensile  failure  was  not  allowed.  This  was 
a preliminary  solution,  intended  primarily  to  exercise  the 
explicit  shear  crack  model  for  the  first  time.  The  results 
have  some  interest  as  a comparison  for  later  cases.  In  addi- 
tion, the  problem  parameters  are  similar  to  those  in  the  small- 
scale  experiments  being  conducted  at  Stanford  Research  Insti- 
tute^. Hence  it  offers  some  degree  of  code  verification. 

The  two  most  prominent  phenomena  which  are  observed  in 
the  SRI  experiments  are  (1)  "flaking"  at  the  springlincs,  and 
(2)  long  vertical  cracks  emanating  from  the  springline  region, 
as  sketched  in  Figure  6. 

Case  H was  run  to  5.24  msec,  which  is  about  2 msec  after 
the  stress  wave  strikes  the  top  of  the  tunnel.  Figures  7,  8, 
and  9 show  the  final  grid  configuration,  velocity  field,  and 
principal  stress  field.  Failed  cells  are  denoted  in  the  grid 
by  line  segments  (aligned  with  the  directions  of  the  shear 
cracks)  at  the  cell  centers.  Extensive  shear  fracture  has 
occurred  at  the  springlines  in  Figure  7,  and  flaking  of  the 
tunnel  wall  would  result  if  the  solution  were  continued.  This 
tendency  is  seen  in  the  effects  of  the  stress  relief  afforded 
by  the  shear  cracking  in  the  velocity  field  (F'igure  8).  At 
5.24  msec,  however,  the  shear  cracking  is  not  propagating  out 
of  the  springline  region  to  form  the  long  vertical  cracks  seen 
in  the  SRI  experiments  (as  sketched  in  Figure  6).  Rather,  the 
stress  plot  in  Figure  9 shows  that  significant  tensile  stresses 
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arc  developing  in  the  regions  where  vertical  cracks  should 
initiate.  These  tensile  stresses  are  of  sufficient  magni- 
tude to  produce  cracking  if  a tensile  failure  model  were  pre- 
sent in  this  solution.  This  suggests  that  the  vertical  cracks 
which  are  observed  experimentally  are  caused  not  by  shear 
failure,  as  had  been  previously  assumed,  but  by  tensile  fail- 
ure. The  stress  relief  allowed  when  the  shear  crack  forms  in 
the  springlinc  region  produces  tension  in  the  surrounding 
material,  and  this  tension  can  cause  fractures  which  propagate 
to  form  the  observed  vertical  cracks. 

4.:  CONCRETE  LINERS  IN  TUFF  - COMPARISONS  SHOWING  EFFECTS 

OF  DIFFERENT  FAILURE  MODELS 

Cases  K and  L examined  the  response  of  concrete  liners 
in  NTS  tuff.  Case  K contained  a 1-ft  linear  elastic  concrete 
liner  in  tuff  with  the  plasticity  failure  model.  Case  L con- 
tained the  same  liner,  but  the  explicit  shear  crack  failure 
model  was  used  for  the  tuff.  Comparison  of  the  results  from 
Cases  K and  L provides  some  indication  of  whether  or  not  use 
of  the  more  refined  explicit  crack  model  is  necessary. 

Figure  10  shows  the  grid  used  for  these  cases.  The  con- 
crete liner  is  denoted  by  the  shaded  region.  Figure  11  shows 
the  computational  grid  in  Case  L at  10  msec,  after  the  tunnel 
has  been  completely  enveloped  by  the  1 kb  stress  wave.  Cells 
which  have  developed  shear  cracks  are  denoted  by  the  crossed 
line  segments  in  the  cells,  with  the  orientation  of  the  line 
segments  indicating  the  directionality  of  the  shear  cracks. 

It  is  seen  that  most  of  the  tuff  cells  in  the  problem  have 
failed. 

Figures  12  and  13  compare  the  stress  fields  for  Cases  K 
and  L at  10  msec.  Both  of  the  failure  models  have  resulted  in 
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Figure  12.  Case  K - Concrete-Lined  Tunnel  in  Tuff 
witli  Plasticity  Failure  Model. 
Principal  Stress  Field  at  10. I msec. 
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Figure  13.  Case  L - Concrete-Lined  Tunnel  in  Tuff 

with  Explicit  Shear  Crack  Failure  Model. 
Principal  Stress  Field  at  10.1  msec. 
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essentially  the  same  stress  fields.  llie  concrete  liner  in 
both  cases  has  been  compressed,  producing  a tangential  (hoop) 
stress  of  about  3 kb.  The  principal  stresses  in  the  tuff  are 
about  1 kb.  Note  that  the  principal  stress  components  in  the 
tuff  adjacent  to  the  liners  are  approximately  equal,  indicat- 
ing that  the  liners  are  being  subjected  to  a uniform  pressure 
around  their  perimeters. 

Figures  14  and  15  show  the  velocity  fields  for  the  two 
lined  cases  at  10  msec,  and  here  we  see  that  the  two  methods 
of  treating  failure  have  produced  significant  differences  in 
the  results.  In  Case  K (Figure  14),  the  tunnel  and  backpack- 
ing have  acquired  a fairly  uniform  velocity,  about  the  same 
as  that  of  the  surrounding  tuff.  In  Case  L (Figure  15),  the 
explicit  shear  crack  model  has  resulted  in  distinctly  non-uni- 
form  velocities.  In  particular,  note  the  large  inward  velo- 
cities on  the  tunnel-liner  interface  at  about  the  10:30  and 
1 : 30  positions . 

These  comparative  results  indicate  that  the  plasticity 
and  explicit  shear  crack  models  for  material  failure  may  pro- 
duce substantially  different  local  velocities  and  hence  dis- 
placements near  the  liner.  In  particular,  the  introduction  of 
explicit  shear  cracks  appears  to  lead  to  non-uniformities  in 
the  velocity  field  near  the  liner.  This  may  cause  stress 
concentrations  which  will  need  to  be  considered  in  the  design 
of  liners.  To  properly  determine  these  concentrations,  mater- 
ial models  and  numerical  methods  which  allow  the  generation 
of  discontinuities  (i.e.,  both  shear  and  tensile  cracks)  will 
be  needed. 

4.3  UNhINED  VS  CONCRETE-LINED  TUNNELS  IN  TUFF 

Case  J consisted  of  an  unlined  tunnel  in  tuff.  The  re- 
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Figure  14.  Case  K - Concrete-Lined  Tunnel  in  Tuff 
with  Plasticity  Failure  Model. 

Velocity  Field  at  10.1  msec. 
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suits  can  be  directly  compared  with  the  preceding  concrete- 
lined  cases  in  tuff  (Cases  K and  L)  to  show  the  action  of  the 
concrete  liner. 


Figures  16  and  17  show  the  stress  and  velocity  fields 
for  the  unlined  Case  J at  the  same  time  as  the  corresponding 
plots  for  the  lined  cases  (Figures  12-15).  The  unlined  tun- 
nel, rather  than  being  accelerated  to  the  velocity  of  the 
surrounding  tuff,  has  developed  significant  inward  velocities 
(at  about  twice  the  velocity  of  the  surrounding  tuff)  and  is 
collapsing.  The  strong  tangential  stresses  which  developed 
in  the  concrete  liner  (e.g..  Figure  12)  are  totally  missing 
in  the  tuff  around  the  unlined  tunnel  (Figure  16) . This  tan- 
gential compression,  of  course,  is  what  prevents  the  collapse 
of  lined  tunnels. 
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Figure  16.  Case  J - Unlincd  Tunnel  in  Tuff 

witli  I’lasticity  Model.  Principal 
Stress  Field  at  10.9  msec. 
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Figure  17.  Case  J - Unlined  Tunnel  in  Tuff 
with  Plasticity  Failure  Model. 
Velocity  Field  at  10.9  msec. 
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5.  COUPLLD  FINITE  DIFFERENCE/FINITE  ELEMENT 
LAP  CODE  FOR  ANALYSIS  OF  THIN  LINERS 

Thick  liners  can  be  represented  using  a finite  difference 
code,  as  is  seen  by  the  examples  in  Cases  K and  L in  Section 
4.2  above.  It  is  not  practical,  however,  to  treat  thin  liners, 
such  as  steel,  by  the  finite  difference  method  because  relative- 
ly small  cells  must  be  used  to  resolve  action  in  the  thin  steel 
section.  The  combination  of  small  cells  and  high  sound  speed- 
in  steel  lead  to  time  steps  which  are  too  small  for  fairly 
long-duration  buried  structure  analyses. 

For  convenience,  steel  liners  are  sometimes  represented 
as  rigid  bodies  in  buried  structure  analyses.  This  was  the 
approach  taken  in  Phase  I of  this  study.  It  is  a reasonable 
approach  so  long  as  clastic  distortion  of  the  liner  does  not 
cause  substantial  stress  adjustments,  or  allow  small  joint 
motions  in  the  nearby  medium.  Unfortunately,  test  results 
indicate  that  elastic  diametral  strains  are  indeed  sufficient 
to  be  of  concern  in  analyses  of  integral  liners.  Thus  blocks 
which  are  in  contact  witli  integral  liners  may  be  able  to 
shift  slightly,  producing  non-uniform  stresses  along  the 
mcdium/liner  interface,  and  perhaps  substantial  stress  con- 
centrations in  the  liner.  These  will  need  to  be  taken  into 
consideration  in  the  design  of  integral  liners. 

Our  approach  to  this  problem  has  been  to  implement  a 
coupled  finite  d i f f erence/ f ini te  element  numerical  method, 
which  we  will  call  LAP,  for  Linked  Analysis  program.  In  this 
coupled  method,  the  earth  medium  is  represented  by  the  Lagran- 
gian  finite  difference  WAVE-L  code,  which  allows  joint  slippage 
and  stress  accommodations  to  occur  around  the  liner,  and 
which  treats  the  far  field  and  near  field  ground  motions  in 
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efficient  manner.  The  integral  structural  liner  is  represent- 
ed by  an  adaptation  of  the  finite  element  SAP  IV  code^,  which 
can  describe  the  displacements  of  the  liner  and  the  dynamics 
of  the  stress  field  in  the  liner.  Concrete  backing  for  a 
steel  liner  can  be  treated  using  either  the  finite  element  of 
the  finite  difference  portion  of  the  LAP  method. 

The  finite  difference  and  finite  element  portions  of 
LAP  are  coupled  on  a cycle-by-cycle  basis  across  a specified 
interface  (usually  the  medium/liner  interface).  In  each 
integration  cycle,  the  finite  difference  code  provides  forces 
across  the  interface  to  the  finite  element  code,  which  in  turn 
provides  displacements  back  across  the  interface.  This  linked 
LAP  method  is  now  operational,  and  it  was  used  to  generate 
the  solutions  with  steel  liners  (Cases  M and  N)  which  are 
described  in  the  following  section. 
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6.  DUl-ORMABLE  STEEL  LINERS 


ft 


For  the  initial  LAR  code  analyses,  two  solutions  were  ob- 
tained for  steel  liners  in  tuff  (Cases  M and  N) . In  both  of 
these  cases,  two- inch  thick  steel  liners  were  used  in  a 10-ft 
dia  tunnel  in  tuff.  The  steel  was  considered  to  be  linearly 
elastic  (i.e.,  with  no  yield  or  failure  criterion).  The  tuff 
was  modeled  using  the  explicit  shear  crack  failure  model  des- 
cribed in  Section  5. 

The  liner  was  represented  in  the  finite  element  portion 
of  the  LAP  code  by  32  beam  finite  elements,  coupled  to  the 
finite  difference  cells  at  the  interface  node  points. 

6.1  TWO- INCH  STEEL  LINER  IN  TUFF  (CASE  M) 

In  tliis  problem,  the  steel  liner  was  in  direct  contact 
with  the  tuff  medium.  The  solution  was  run  to  15  msec.  Fig- 
ures 18,  19,  and  20  show  the  computational  grid  velocity  field, 
and  principal  stress  fields  at  7.5  msec,  as  the  stress  wave  is 
engulfing  the  tunnel.  Figures  21,  22,  and  23  show  the  same 
plots  at  14.5  msec,  with  the  tunnel  completely  engulfed  by  the 
step  wave.  The  pattern  of  shear  cracks  is  similar  to  that  ob- 
tained in  the  earlier  solution  of  a concrete  liner  in  tuff 
(Case  L , F igure  11). 

Time  histories  of  the  steel  liner  response  were  obtained 
from  the  finite  element  portion  of  LAP  at  the  stations  identi- 
fied in  Figure  24.  Figures  25  to  29  show  the  hoop  stress 
histories  on  the  inner  surface  (i.e.,  the  free  surface)  and 
on  the  outer  surface  (the  steel/tuff  interface).  These  stresses 
exceed  a nominal  yield  strength  of  steel  (say  10  kb,  or  150  ksi) 
after  about  12  msec,  so  permanent  deformation  and  probably 
failure  of  the  liner  will  occur.  At  the  end  of  the  solution 
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at  15  msec,  the  hoop  stress  is  beginning  to  level  oft'  as  it 
approaches  about  30  kb  (450  ksi)  , which  would  be  the  stress 
level  reached  in  a static  solution  of  this  problem. 

Figures  30  to  35  show  displacement  and  velocity  histories 
at  the  top  and  bottom  stations,  and  at  the  springline  of  the 
liner.  The  vertical  displacements  and  velocities  are  relative 
to  the  center  of  mass  of  the  steel  liner.  Note  that  the  top 
station  (25)  has  been  pushed  down  towards  the  center  of  mass 
about  4 inches  by  15  msec.  The  bottom  station  (9)  first  moves 
up,  relative  to  the  mass  center,  then  also  moves  downward,  to 
about  1 inch  by  15  msec.  The  springline  station  (1)  exper- 
iences about  1-in  outward  radial  displacement  by  15  msec  (about 
2- in  diametral  strain).  These  motions  occur  at  velocities, 
relative  to  the  center  of  mass,  in  the  range  up  to  about  150 
fps.  Slopes  of  the  velocity  histories  indicate  relative 
accelerations  of  the  order  of  1000-3000  g's. 

6.2  TWO- INCH  STEEL  LINER  BACKED  BY  TWO  FEET  OF  CONCRETE  IN 

TUFF  (CASE  N) 

This  problem  was  identical  to  the  preceding  Case  M,  except 
that  the  steel  liner  was  backed  by  about  two  feet  of  concrete. 
The  concrete  was  treated  by  the  finite  difference  portion  of 
the  code,  so  the  finite  difference/finite  element  interface 
was  between  the  steel  liner  and  the  concrete.  Figure  36  shows 
the  computational  grid  at  7.5  msec,  as  the  stress  wave  was  en- 
gulfing the  tunnel.  The  concrete  backing  is  shown  by  shading 
in  this  plot.  Figure  37  and  38  snow  the  velocity  and  princi- 
pal stress  fields  at  this  same  time. 

The  solution  was  carried  out  until  9.5  msec,  and  the  final 
grid,  velocity  field,  and  principal  stress  field  are  shown  in 
Figures  39,  40,  and  41,  The  stress  level  in  the  concrete  had 
reached  about  3-4  kb  by  this  time,  and  was  beginning  to  rise 
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sharply  in  the  steel.  We  did  not  continue  the  solution  further. 

An  alternative  to  placing  the  finite  di f f erence/ f init e 
element  interface  between  the  steel  and  concrete  would  be  to 
include  the  concrete  in  the  finite  element  portion  of  the 
solution.  This  would  be  somewhat  more  efficient,  but  it  would 
sacrifice  the  generality  in  treatment  of  material  properties, 
fracture,  and  explicit  crack  propagation  which  is  available 
in  the  finite  difference  code. 
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Figure  18.  Cdmputational  Grid  for  Lagranglan  Finite 
Difference  Portion  of  LAP  Code  Solution 
for  Case  M - 2"  Steel  Liner  in  Tuff  with 
l■.xl)lll•lt  SluMr  Crack  Failure  Model. 
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Figure  19.  Velocity  Field  for  Lagranglan  Finite  Difference 

Portion  of  LAP  Code  Solution  for  Case  M - 2"  Steel 
Liner  in  Tuff  with  Explicit  Shear  Crack  Failure 
Model  (t  » 7.5  msec) 
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Figure  20.  Principal  Stress  Field  for  Lagranglan  Finite 
Difference  Portion  of  LAP  Code  Solution  for 
Case  M - 2"  Steel  Liner  in  Tuff  with  Explicit 
Shear  Crack  Failure  Model.  (t  ■ 7.5  csec) 
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Figure  21.  Computational  Grid  for  Lagrangian  Finite 

Difference  Portion  of  LAP  Code  Solution  for 
Case  M - 2"  Steel  Liner  in  Tuff  with  Explicit 
Shear  Crack  Failure  Model  (t  » 14.5  msec) 
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Figure  22.  Velocity  Field  for  Lagrangian  Finite  Difference 

Portion  of  LAP  Code  Solution  for  Case  M - 2"  Steel 
Liner  in  Tuff  with  Explicit  Shear  Crack  Failure 
Model  (t  “ 14.5  msec) 
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Figure  23.  Principal  Stress  Field  for  Lagranglan  Finite 
Difference  Portion  of  LAP  Code  Solution  for 
Case  M - 2"  Steel  Liner  in  Tuff  with  Explicit 
Shear  Crack  Failure  Model.  (t  » 14.5  msec) 
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Figure  2U.  Locations  of  Tinie-Mlstory  Stations  for 

Plots  In  Figures  25-29  (Case  M - 2"  Steel 
Liner  in  Tuff) 
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Figure  25.  Compressive  Hoop  Stress  History  at  Station  25 

in  2"  Steel  Liner.  Case  M - Steel  Liner  In  Tuff 
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Figure  26.  Compressive  Hoop  Stress  History  at  Station  29 

in  2"  Steel  Liner.  Case  M - Steel  Liner  in  Tuff 
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Figure  27.  Compressive  Hoop  Stress  History  at  Station  1 

in  2"  Steel  Liner.  Case  M - Steel  Liner  in  Tuff 
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Figure  28.  Compressive  Hoop  Srress  History  at  Station  5 

in  2"  Steel  Liner.  Case  M - Steel  Liner  In  Tuff 
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Figure  29.  Compressive  Hoop  Stress  History  at  Station  9 

in  2"  Steel  Liner.  Case  M - Steel  Liner  in  Tuff 
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Figure  30.  Vertical  Displacement  (Relative  to  Center  of 
Mass  of  Steel  Liner)  of  Station  25. 

Case  M - 2"  Steel  Liner  in  Tuff 
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Figure  31.  Vertical  Displacement  (Relative  to  Center  of 
Mass  of  Steel  Liner)  of  Station  9. 

Case  M - 2"  Steel  Liner  in  Tuff 
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Figure  32.  Horizontal  (Radial)  Displacement  of  Station  1 

in  2"  Steel  Liner.  Case  M - Steel  Liner  in  Tuff 
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Figure  33.  Vertical  Velocity  of  Station  25  (Relative  to 
Center  of  Mass  of  Steel  Liner).  Case  M - 
2"  Steel  Liner  in  Till  f 
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Figure  3A.  Horizontal  (Radial)  Velocity  of  Station  1 

in  2"  Steel  Liner.  Case  M - Steel  Liner  in  Tuff 
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Figure  35.  Vertical  Velocity  of  Station  9 (Relative  to 
Center  of  Mass  of  Steel  Liner).  Case  M - 
2"  Steel  Liner  in  Tuff 
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Figure  36.  Computational  Crld  for  Lagranglan  Finite 
Difference  Portion  of  LAP  Code  Solution 
for  Case  N - 2"  Steel  Liner  with  2*  Concrete 
Backing  in  Tuff  with  Explicit  Shear  Crack 
Failure  Model  (t  = 7.5  msec) 
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Figure  37.  Velocity  Field  for  Lagrangian  Finite  Difference 
Portion  of  LAP  Code  Solution  for  Case  N - 
2"  Steel  Liner  with  2'  Concrete  Backing  In  Tuff 
with  Explicit  Shear  Crack  Failure  Model 
(t  » 7.5  msec) 
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Figure  38.  Principal  Stress  Field  for  Lagrangian  Finite 
Difference  Portion  of  LAP  Code  Solution  for 
Case  M - 2"  Steel  Liner  in  Tuff  with  Explicit 
Shear  Crack  Failure  Model.  (t  « 7.5  msec) 
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Figure  39.  Computational  Grid  for  Lagranglan  Finite 
Difference  Portion  of  LAP  Code  Solution 
for  Case  N - 2"  Steel  Liner  with  2'  Concrete 
Backing  In  Tuff  with  Explicit  Shear  Crack 
Failure  Model  (t  • 9.5  msec) 
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Figure  40.  Velocity  Fielil  for  Lagranglan  Finite  Difference 

Portion  of  LAP  Code  Solution  for  Case  N - 2"  Steel 
Liner  wltli  2'  Concrete  Racking  In  Tuff  with  Explicit 
Shear  Crack  Failure  Model.  (t  = 9.5  msec) 
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Figure  41.  Principal  Stress  Field  for  Lagrangian  Finite 
Difference  Portion  of  LAP  Code  Solution  for 
Case  M - 2"  Steel  Liner  in  Tuff  with  Explicit 
Shear  Crack  Failure  Model  (t  • 9.5  msec) 
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ATTN:  Code  04B 
ATTN:  Technical  Library 

Superintendent  (Code  1424) 

Naval  Postgraduate  School 

ATTN:  Code  2124,  Tech.  Rpts.  Librarian 
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Commander 

Naval  Surface  Weapons  Center 
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ATTN:  Code  EM,  Jack  Kalinowski 

President 
Naval  War  College 
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ATTN;  DES-S,  M.  A.  Plamondon 
ATTN:  SUL 

Headquarters 

Air  Force  Systems  Conmand 
ATTN:  DLCAW 

Commander 

Foreign  Technology  Division,  AFSC 
ATTN:  NICD,  Library 

HQ  USAF/IN 

ATTN;  INATA 

HQ  USAF/PR 

ATTN:  PRE 

HQ  USAF/RD 
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